function dx=Plant(t,x,flag,para)
dx=0;

id=para(1);
uq=para(2);
w=para(3);

Rs=1.74;np=4;Lq=10;fi=0.353;

dx(1)=-np*w*id-Rs*x(1)/Lq-np*fi*w/Lq+uq/Lq;

